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ABSTRACT 



High resolution spectra of 123 red giant stars in the globular cluster M13 and 64 
red giant stars in M92 were obtained with Hectochelle at the MMT telescope. Emission 
and line asymmetries in Ha, and Ca II K are identified, characterizing motions in 
the extended atmospheres and seeking differences attributable to metallicity in these 
clusters and M15. On the red giant branch, emission in Ha generally appears in stars 
with T e ff < 4500 K and log L/Lq > 2.75. Fainter stars showing emission are asymptotic 
giant branch (AGB) stars or perhaps binary stars. The line-bisector for Ha reveals the 
onset of chromospheric expansion in stars more luminous than log (L/Lq) ~ 2.5 in all 
clusters, and this outflow velocity increases with stellar luminosity. However, the coolest 
giants in the metal-rich M13 show greatly reduced outflow in Ha most probably due to 
decreased T e ff and changing atmospheric structure. The Ca II K3 outflow velocities 
are larger than shown by Ha at the same luminosity and signal accelerating outflows in 
the chromospheres. Stars clearly on the AGB show faster chromospheric outflows in Ha 
than RGB objects. While the Ha velocities on the RGB are similar for all metallicities, 
the AGB stars in the metal-poor M15 and M92 have higher outflow velocities than in the 
metal-rich M13. Comparison of these chromospheric line profiles in the paired metal- 
poor clusters, M15 and M92 shows remarkable similarities in the presence of emission 
and dynamical signatures, and does not reveal a source of the 'second-parameter' effect. 

Subject headings: stars: chromospheres - stars: mass loss - stars: AGB and post-AGB 
- globular clusters: general - globular clusters: individual (M13, M15, M92, NGC 2808) 



Introduction 



The well— known second-parameter problem in globular clusters (jSandage fe Wildevl 119671 ). 
in which a parameter other than metallicity affects the morphology of the horizontal branch, re- 
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mains unresolved. Metallicity, as first noted by ISandage fc Wallersteinl (|1960l ). remains the prin- 
cipal parameter, but pairs of clusters, with the same metallicity, display quite different horizontal 
branch morphologies thus challenging the canonical models of stellar evolution and leading to the 
need for a 'second-pa r ameter 1 . Cluster ages have been examin e d in many studi e s (ISearle &: Z inn 



1978; iLee et all 1 1994 IStetson et al.l Il996l ; iLee fc Carnevi Il999 : ISaraiedinil Il997l ; ISaraiedini et al 



19971 ) and in addition, many other suggestions for the 'second-parameter (s)' have been proposed, 



including: total cluster mass; stellar environment (and possibly free-floating planets); primor- 



mass loss ( Catelan 


2000: 


Catelan et al. 


2001 


; Sills & Pinsonneaultl 


2000: 


Soker ct al. 


2001; 


Sweieartl 


1997; Buonanno et al. 


1993: 


Peterson et al. 


1995: 


Buonanno et al. 


1998: 


Recio-Blanco et al.l 12006). 


Many authors ( 


VandenBerg et al. 


1990: 


Lee et al. 


1994: 


Catelan 


2000) have proposed that the 



second-parameter problem cannot be explained by only one second parameter. Various studies 
have attempted to limit the number of second parameters. 

An example of paired second-parameter clusters is M15 and M92 ([Fe/ H|=— 2.26 and —2.2 8 



respectively). Although the metalliciti es of these two cluster s are the same (ISneden et al.ll2000l ). 
their horizontal branches (HBs) differ (jBuonanno et al.lll985l ). M92 has a brighter (by about one 
magnitude) and redder blue HB exte nsion than Ml 5. T he color magnitude diagrams (CMDs) of 
this pair were examined in detail by ICho &: Led (|2007l ). They found that the difference in the 
HB morphology between the two is probably not a result of deep mixing in their red giant branch 
(RGB) sequences, because no significant 'extra stars' were found i n their observed RGB luminosity 
functions compared to the theoretical RGB luminosity functions. ISneden et al.l (|2000l ) found that 
Si, Ca, Ti, and Na abundance ratios of the red giants are nearly the same in both clusters, only the 
[Ba/Ca] ratio shows a large scatter and the mean value in M15 is twice that found in M92. These 
studies eliminate deep mixing and subtle abundance variations as possible second parameters, while 
m ass loss is examined in this paper. Detailed observations of red giant stars in M15 are contained 
in 



Meszaros et al.l (|2008l ). but the comparison between M15 and M92 is described here. 



M13 ([Fe/H]=— 1.54) is one of the most studied second-parameter globular clusters. M13 and 
M3 are almost identical in most respects (metallicity, age, chemical composition), but there are dra- 
matic differences in both the HB and blue straggler populations. Analysis of both clusters' CMDs 
(jFerraro et al.lll997l ) with the Hubble Space Telescope revealed that neither age nor cluster density, 
two popular second-parameter candidates, is likely to be responsible for the differences in these 
clusters. From the analy sis of high-resolut i on, h igh signal-to-noise ratio spectra of six RGB stars 
in M3 and three in M13, ICavallo Nagarl ()200ol ) found that the [Al/Fe] and [Na/Fe] abundances 
increase toward the tip of the RGB. They conclu ded that the da t a for both clusters are consistent 
with deep mixing as a second parameter. Later, iJohnson et al.l (|2005l ). from medium— resolution 
spectra of more than 200 stars in M3 and M 13, concurred that deep m ixing is the best candidate 
for second parameter in this pair of clusters. ICaloi D'Antonal (|2005l ) also examined the second- 
parameter problem in M3 and M13 in detail and proposed that the overall difference between M3 
and M13 CMD morphologies is due to the different helium content. Since M13 does not have a red 
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clump in its horizontal branch they suggested that it represents an extreme case of self-enrichment 
of helium, which might come from the massive asymptotic giant branch stars (AGB) in the first 
~100 Myr of the cluster life. 



A multivariate study of the CMDs for 54 globular clusters was carried out by lRecio-Blanco et al 



([200a) from Hubble Space Telescope photometry to quantify the parameter dependencies of HB 
morphology. They found that the total cluster luminosity ( theref ore the total mass) has the largest 
impact on the HB morphology, and as ICaloi D'Antonal (120051) speculate d, there may be enrich- 
ment of helium from an earlier population of stars. iD'Antona et al.l (|2002l ) modeled the evolution 
of globular cluster stars and showed that different choices of mass— loss rate affect the distribution 
of stars on the HB. 

In this paper we discuss the characteristics of Ha and Ca II K emission in M13 ([Fe/H]=— 1.54) 
and M92 ([Fe/H]=— 2. 28). We compare our results with the previously observed, metal-poor cluster 
M15 ([Fe/Hl=-2.26) (jMeszaros et al]l2008l ) and the metal-rich cluster NGC 2808 ([Fe/H]=-1.15) 



(jCacciari et al.ll2004l ). Detailed study of these four clusters allows us to examine a possible depen- 
dence between the average cluster metallicity and characteristics of Ha and Ca II K emission, and 
diagnostics of mass outflow. Observations with the same instrument of the second-parameter pair 
M15 and M92 offer a good comparison to examine mass loss as a possible second parameter. 



Observations and Data Reduction 



The Hectochelle on the MMT (jSzentgyorgyi et al.lll998l ) contains 240 fibers that can be placed 
~ 2 arcsec apart on the sky across the field of view which spans 1 degree. The diameter of one 
fiber on the sky is 1.6 arcsec. The apparent diameter of M13 is ~15 arc minutes, and about 60—70 
red giants in the globular cluster could be measured with each fiber configuration. The apparent 
diameter of M92 is smaller, and only 30—40 stars could be measured with one configuration. Two 
separate input fiber configurations for different stars were made for each cluster. A total of 123 
different red giant stars in M13 and 64 red giants in M92 were observed in 2006 March and 2006 
May. 

Targe t s wit h a high probabilit y (> 95%) of membership were chosen from the catalog of 
Cudworthl (|1976T ) for M13 and from ICudworth fc Monetl (119791 ) for M92. Smooth coverage of the 



RGB and AGB could be achieved within the constraint of the fiber placement on the sky. The 
CMD of the observed cluster members can be seen in Figure 1 and the target stars are listed in 
Table 1 (M13) and Ta ble 2 (M92). Coordinates of the stars were taken from the 2MASS catalog 
(jSkrutskie et al.1 120061 ) and used to position the fibers. Many fibers (~ 150—200) were placed 
on blank regions of the sky in order to measure the sky background in detail. These sky fibers 
were equally distributed in the Hectochelle field of view to cover a large area around the clusters. 
Hectochelle is a single-order instrument and three orders were selected for observation with order- 
separating filters: OB25 (Ha, region used for analysis AA 6475 — 6630), Ca41 (Ca II H&K, region 
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used for analysis AA 3910 — 3990), and RV31 (region used for analysis AA 5150 — 5300). Bias and 
quartz lamps for the flat correction were obtained during the afternoon each day. Exposures with 
the ThAr comparison lamp were obtained before and after every observation during the night to 
determine the wavelength solution. The spectral resolution was about 34,000 as measured by the 
FWHM of the ThAr emission lines in the comparison lamp. Exposures in each filter are summarized 
in Table 3. The number of objects observed changed slightly between observations, because fiber 
positions need to be reconfigured when targets pass the meridian. 

Data reduction was done using standard IRAF spectroscopic packages. The uneven sky inten- 
sity in the CCD required special n on-standard methods of sky subtraction. These procedures are 
described in lMeszaros et al.l f)2008l ) . 



3. Stars with Ha Emission 

3.1. Determining Ha Emission 

Two methods have been used for the identification of Ha emission. The first selects stars 
with flux in the Ha wings lying above the local continuum. Strong emission can be easily found, 
however faint emis sion comparable to th e noise of the continuum can be missed. A second method 



was introduced bv lCacciari et al.l (|2004l ) and is illustrated in Figure 2. They select a star without 
emission and the Ha absorption line from this star is subtracted from the other spectra. With this 
method weak emission can be identified, but it strongly depends on the template selected. The 
Ha absorption profile depends on temperature, as well as broadening from turbulent velocity and 
rotation, both of which could introduce features in the subtracted profile. An individual Kurucz 
model can be made for every star as a template, and the temperature problem can be avoided, 
but the uncertainty of other physical parameters can introduce similar effects in the subtracted 
spectrum. In a continuum— normalized spectrum, the Ha emission appears above the continuum 
level for the majority of stars (L465 and L72 in Figure 2). However in fainter stars only the shape 
of the Ha line profile is changing and the emission does not appear above the continuum, rather 
just a small additional flux emerges in the absorption wings (L250 and L252 in Figure 2). The 
identification of this kind of emission can be challenging (L403 in Figure 2). In this paper, we 
used both methods. However, we selected 8 stars with no emission and of different colors and 
luminosities to make the template. The stars identified with Ha emission are the same with both 
methods. 

No matter which method is used, the detection of faint emission depends on the reduction 
technique. Continuum normalization and sky subtraction can change the emission flux and move it 
above or below the continuum level (L252 in Figure 2). Continuum normalization was done using 
a low— order Chebyshev function in the IRAF task continuum in order to fit the continuum and 
filter throughput. The continuum placement strongly depends on the order of the function and the 
rejection limits below and above the fit. Sky subtraction is especially challenging with Hectochelle 
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because additional counts appear between aperture numbers 100 and 150, possibly due to scattered 
light. This additional flux depen ds not only on the ap erture but also wavelength, and although a 
reduction system was developed (jMeszaros et al.ll2008l ). all sky background cannot be subtracted 
in the middle section of the CCD. For these reasons the emission of very faint stars can be hard to 
identify and this can introduce uncertainties in the statistics of the presence of emission. 



3.2. Ha Emission in the CMD 



Emission in Ha signals an extended and high— temperatur e chromosphere; in addition the 



asym metry of the emission indicates chromospheric mass motions (jDupree et al.lll984l ; iMauas et al 



2006). We observed a total of 123 different red giant stars in M13 and found 19 with Ha emission. 
In M92, we found 9 stars with Ha emission out of 64 objects. Emission above the continuum in the 
Ha profile can be seen in Figures 3 and 4. For comparison, Figure 4 includes a star that exhibits no 
emission. The color-magnitude diagram (CMD) for each night of observation appears in Figure 5 
for both clusters. The intensity ratio, B/R, of Blue (short wavelength) and Red (long wavelength) 
emission peaks for stars showing emission is contained in Table 4. 

In M13, emission is found in stars brighter than V= 1 4.69, correspo nding to My = +0.21, using 
the apparent distance modulus (m — M)y = 14.48 from lHarrisI (|1996l ). The star which marks the 
faint luminosity limit (L719) appears to be either a blend, or a physical binary RGB star judging 
from its position in the CMD (bluer than RGB stars at the same absolute magnitude). This star 
also had a significant radial velocity change between observations (see Section 4.1). Among the 
RGB stars, L1073 at V=12.88 {My = —1.60) marks the faint luminosity limit of Ha emission. 
Stars brighter than this are on the RGB or AGB; the evolutionary status cannot be determined 
from the CMD itself. Stars in M92 that show emission are br ighter t han V =14.54 {My = —0.1), 
using the apparent distance modulus {m — M)y = 14.64 from lHarrisI (1996). However the faintest 
star (IX-12) showing emission appears to be an AGB star, according to its position in the CMD. 
Considering stars on the RGB, the star IV-94 (V=13.06, My = —1.58) appears to be the faintest 
RGB star showing emission in M92 (although the differences between the RGB and AGB at that 
part of the CMD are very small). In M15, the faint luminosity limit showed significant changes 
between observations ; this amounted to a change in the faint magnitude limit of 0.79 magnitudes 
(jMeszaros et al.ll2008l ). One can assume that the emission behaves very similarly in these clusters 
as well, and that the faint luminosity limit of Ha emission is not constant. 

For comparison on a luminosity scale, unreddened colors for M13 and M92 stars were cal- 
culated. Foreground red dening \ E(B — V) = 0.02 for both clusters] and the apparent distance 
modulus were taken from lHarrisI (|1996l ). The effective temperatures, bolometric corrections, and 
luminosities were obtaine d from the V — K colors (Tables 5 and 6) using the empirical calibrations by 
Alonso et all dl999l. boOlh and the cluster average metallicity [Fe/H]=-1.54 for M13, [Fe/H]=-2.28 
for M92 (jHarriall996l ). Thus on the red giant branch alone, emission appears in stars brighter than 
log {L/Lq) = 2.79 in M13 and ~ 78% of these stars (18) show Ha emission. In M92 this luminosity 
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limit is slightly lower than in M13, log (L/Lq) = 2.74, and also ~ 78% of these stars (7) show 
emission. 

Although both clusters were observed on two different days, the configurations were chosen to 
eliminate stars already observed in order to achieve full coverage of stars in the CMD. When it was 
possible, previously observed stars were configured to the remaining fibers, but the number of stars 
observed twice for both clusters is very small, only 17 in M13 and 15 in M92. Of the stars showing 
Ha emission, comparison was possible for only two stars in M13 and three stars in M92. In M13, 
between 2006 March 14 and 2006 May 10, L72 changed asymmetry (see Figure 3), while for the 
other star, L719, the already weak emission vanished. In M92, all three stars (11-53, VII-18, and 
IX-12) kept the same emission asymmetry, but the flux level of IX-12 changed in only two days 
(see Figure 4). 



4. Radial Velocities 

4.1. Cross-correlation Technique 

To measure accurate radial velocities, we chose the cross-correlation method using the IRAF 
task xcsao. Two filter regions, OB25 and RV31, were used for radial velocity measurements. The 
spectral region on the RV31 filter between 5150 A and 5300 A contains several hundred narrow 
photospheric absorption lines of predominantly neutral atoms and very few terrestrial lines, thus 
the cross-correlation function is narrower than from the Ha region, which only contains ~10 lines 
(Figure 6). In the OB25 filter, the region selected for the cross-correlation spanned 6480 A to 
6545 A purposely omitting the Ha line. This results in 100—200 m s _1 errors with the RV31 filter as 
compared to 200—400 m s -1 using the wavelength region earlier described in the Ha filter. Spectra 



of our targets from both filters were cross-correlated against 2280 spectra calculated bv lCoelho et al 



(|2005l ) covering temperatures between 3500 and 7000 K, metallicities between [Fe/H]=— 2.5 and 
+0.5, and log g between and 5. Radial velocities were corrected to the solar system barycenter. 
To calculate the radial velocity of a star, the radial velocities from ten templates with the highest 
amplitude of the cross-correlation function for each filter were collected and averaged together. A 
sample of the template spectra compared to an observation can be seen in Figure 6. The physical 
parameters of the templates that were used for the radial velocity measurements usually agreed 
with each other within 200 K in temperature, 1 in log g, and —0.5 in [Fe/H] with our calculated 
physical parameters (Tables 5 and 6). For almost every star the radial velocity differences among 
the 10 highest correlation templates in each filter were less than 0.5 ± 0.2 km s~ x , which is close 
to the error of the individual measurements. 



We compare our results with those found in the literature. In M13, ISoderberg et al.l (| 19991 ) 
used the Hydra spectrograph on the 4— m Mayall telescope to obtain spectra of 150 stars. Their 
template for the cross-correlation was an averaged spectrum of all giants for each Hydra observa- 
tion. Therefore the individual radial velocities were determined as compared to the average cluster 
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velocity. The radial velocity of the averaged spectrum was calculated by cross-correlating it to the 
solar spectrum. Comparison of the results can be seen in Figure 7. Errors spanned 0.5 km s _1 to 
3—5 km s~ x in their sample, and there is a system atic 1.1 ±0.5 km s -1 of fset (Figure 7, left upper 
panel) between our radial velocities and those of ISoderberg et al.l (|1999l ). Hectochelle velocities 
determined using the Ha region from 2006 March 14 agreed with the observations two days later 
with the RV31 filter (see Figure 7, left lower panel) for the same stars. Radial velocities calculated 
from the data taken with the RV31 filter in 2006 May also agreed with data taken with the OB25 
filter on 2006 March 14 (Figure 8, right upper and lower panel). We find the average radial ve- 
locity of M13 to be —243.5 ± 0.2 k m s -1 , which i s slightly lower than the cluster radial velocity 
(—245.6 ±0.3 km s" 1 ) quoted in the lHarrisI (|1996l ) catalog. 



Five stars observed with Hectochelle in M13 were reported as possible binaries by IShetrone 



(| 19941 ). when the radial velocities measure d with the 3-m Shan e telescope (Lick Observatory) were 
compared with velocities determined by iLupton et al.l (jl987l ). In all of these stars, differences 
between the two observations were larger than 4 km s _1 which exceeds the measurement errors 
of ~ 1 km s _1 and may reflect intrinsic stell ar variability o r bin ary reflex motions. Our radial 
veloc ities differ by 4 — 5 km s _1 compared with lLupton et al.l ( 19871 ). but agree within 1 — 2 km s _1 
with iShetrond (| 19941 ). which also suggests that long-term changes are present. Among these five 
stars, we o bserved one, L72, w hich showed 2.1 km s _1 velocity change between 2006 March and 
2006 May. ILupton et al.l (|1987l ) identified this star in M13 as a possible binary from variations in 
radial velocity over several years of observations. L72 is al so known as a pulsating variable star with 
a possible period of 41.25 days (IRusseva fc Russevlll980l ). so the velocity change found here may 
also relate to pulsation. L719, which marks the faint luminosity limit of stars showing Ha emission, 
also had radial velocity changes between 2006 March 14 to 2006 May 10 from —254.1 ±0.3 km s _1 
to -245.2 ±0.2 km s -1 . 
sample in M13. 



No other stars showed significant, larger than 2 km s , variations in our 



A large sample of stars in M92 was observed bv lDrukier et al.l (|2007l ) using the HYDRA multi- 
fiber spectrograph on the 3.5-m WIYN telescope. Their errors spanned 0.3—1.2 km s . The 
comparison of results can be seen in Figure 8. Radial velocities for the same stars agreed within 
the errors (Figure 8, left upper panel). The Hectochelle spectra give the cluster average radial 
velocit y as — 118.0 ±0.2 km s , which is lower than the value (—120.3 ±0.1 km s -1 ) quoted in the 
Harris! (119961 ) catalog. In M92, two stars show radial velocity variations, which usually indicates 
binarity or pulsation. 11-53 had a significant velocity variation of 7.7 km s _1 between 2006 May 7 
and 2006 M ay 9. Another s tar, X I-38, showed a 4.9 km s _1 difference between the radial velocity 
measured bv lDrukier et al.l (|2007l ) and the velocity measured by us on 2006 May 7. 



x If this object were a single line binary, the velocity change allows only lower limits to the period (P>90 days) 
and semi-major axis (a >10R*). The putative companion to the red giant could be either a white dwarf or a late 
main sequence star, and probably the former since the color is bluer than a red giant. 
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4.2. Bisector of Ha Lines 

To search for mass motions in the chromosphere, we evaluated the Ha core asymmetry using 
a bisector method. The difference between the center of the line core and the line center near 
the continuum level gives a measure of the atmospheric dynamics through the chromosphere. To 
accomplish this, the line profile with continuum normalization was divided into 20 sectors. The top 
sector was usually between 0.7 and 1.0 of the continuum in the normalized spectrum, the lowest 
sector was placed 0.01 — 0.05 above the deepest part of the line depending on its signal-to-noise 
ratio. The velocities of the Ha bisector asymmetry (vbis) are calculated in the following way: the 
top and the bottom 3 sectors are selected, the wavelength average of each sector is calculated, then 
subtracted one from another and changed to a velocity scale. The bisector velocities, v^is, are shown 
in Figure 9 and listed in Tables 7 and 8. A negative value corresponds to an outflowing velocity. The 
error of the majority of the measurements spanned 0.5—1.0 km s _1 ; only stars fainter than V=15 
magnitude exceeded 1.0 km s _1 in measurement error. RGB stars fainter than log (L/Lq) = 2.5 
did not show asymmetry in the Ha core and vyis was nearly zero. This lu minosity is nearly the 



same for both clusters and also very similar to M15 (IMeszaros et al.ll2008l ). which suggests that 
the luminosity limit of the line core asymmetry marking the onset of expansion does not depend 
on average cluster metallicity. Stars brighter than log (L/Lq) = 2.5 showed core asymmetry and 
the majority of the bisectors were blue shifted. The start of chromospheric outflow presumably 
relates to mass loss. However, the value of vus appears to depend on luminosity. In the metal-rich 
cluster M13, v^s increases with luminosity and reaches its maximum value (~ 5 km s _1 ) at about 
log (L/Lq) = 3.2 but the most luminous stars exhibit lower (near zero km s _1 ) values. In the 
more metal-poor cluster M92, v^s also increases with luminosity and reaches the maximum value 
at log (L/Lq) = 3.4 but the decrease in outflow velocity is much smaller. At the same luminosity, 
T^eff f° r the metal-rich M13 giants is lower than for M92. Thus, the apparent decrease in at 
high luminosity would appear to be related to the changing atmospheric structure (see Section 6.1). 

Where AGB stars are well separated in color from the RGB stars in M92, the AGB stars 
exhibit larger values of v^s than RGB stars (Figure 10). The star IX-12 in M92 shows the largest 
bisector velocity at Vbis =—15.9 ± 1.3 km s" 1 , and this star is the faintest star showing emission in 
the cluster (see Figure 9). Its position in the CMD suggests that this star is an AGB star (Figure 5). 
AGB stars in M13 generally have lower bisector velocities than AGB stars in M92 which suggests 
that the metal-rich objects have slower winds. 



5. Ca II K Profiles 

Spectra in the Ca II H&K region were obtained for 119 red giant stars in M13 and 63 red 
giants in M92. The profiles of the Ca II K core (A3933) are shown in Figures 11 and 12 for all stars 
exhibiting emission. The position of these stars in the CMD can be seen in Figure 13. The intensity 
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ratio of the emission corej^l B/R, is summarized in Table 9. Because of the high radial velocity of 
M13 {v ra d=— 243.5 km s _1 ), absorption by the local interstellar medium (ISM) is well away from 
the Ca II K core. Although the ISM is closer to the Ca II K emission in M92 (v ra d=— 118.0 km s _1 ), 
it does not affect the emission profile. Red giant stars have low flux levels near 3950 A, hence the 
deep photospheric absorption in Ca II H&K causes the photon noise to become comparable to 
the flux of the core emission for stars fainter than 14th magnitude. Determining the presence of 
e mission is challenging . The spectra of our targets were compared to a Kurucz model, [computed 



by lCoelho et al.1 (|2005l ) without a chromosphere], to verify the emission. Altogether 34 stars showed 



Ca II K emission in M13 and 12 in M92. 

The spectra of Ha and Ca II were obtained on the same night or separated by 1 or 2 days, 
and the asymmetry of the K-line emissions is similar to that found in Ha for the majority of the 
stars. The brightest stars in M13 showed stronger Ca II K emission than stars in M92 at the 
same luminosity. This results because the lower effective temperatures of M13 giants increase the 
contrast of the emission to the continuum. For stars fainter than V=14, the ratio of the blue to the 
red side of the Ca II K emission core, or even the presence of the emission is difficult to determine. 

The velocity of the central reversal (K3) was measured for the brightest stars (Table 10) using 
three strong absorption lines closest to Ca II K as a photospheric reference. A Gaussian function 
using the IRAF task splot was fitted to the cores of the photospheric lines and the central reversal 
absorption of Ca II K. Radial velocities of the photospheric lines were averaged and then subtracted 
from the radial velocity of the Ca II K3 feature. The velocity shift of the K3 absorption lies between 
km s" 1 and —16 km s -1 (Table 10). 



6. Discussion 



In this paper and in our previous study (jMeszaros et al.l 120081 ) . we have prese nted Ha an d 
Ca II K spectroscopy of 297 red giants in 3 globular clusters with different metallicities ( Harrisfl996 ): 
M13 ([Fe/H]=-1.54), M15 ([Fe/H]=-2.26), and M92 ([Fe/H]=-2.28). The presence of emission 
in these transitions signals an extended, high-temperature chromosphere, and the asymmetry of 
the emission and the line core indicates chromospheric mass motions. Comparison of the statistics 
of the profile characteristics among the globular cluster stars could reveal the effects of metallicity 
on mass loss. In the following sections, we discuss several characteristics of the presence of Ha and 
Ca II emission and the resulting velocity signatures. Parameters of the clusters and the lines can 
be found in Table 11 . In the final section of the Discussion, we compare our results with those of 
Cacciari et al.l (|2004l ) who presented similar line profiles for 137 red giants in the globular cluster 
NGC 2808. 



2 B signifies the short-wavelength emission peak and R the long-wavelength emission peak. 
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6.1. Presence of Ha Emission 

On the RGB, Ha emission sets in for all stars with T e ff < 4500 K and log (L/Lq) > 2.75 in 
all 3 clusters: M13, M15@ and M92. It is perhaps fortuitous that the limits are so similar since 
the prese nce of Ha can chang e by as much as 0.79 magnitudes from observations on one date to 
another ( Meszaros et alJbood ). Stars on the AGB show Ha emission to lower luminosity limits 
than the RGB objects. The faint limits of emission for AGB stars in M13 and M92 are comparable, 
while AGB stars in M15 with emission are brighter. Emission is variable in all giants and again 
there does not appear to be a systematic dependence of luminosity limits on metallicity. This result 
suggests that whatever mechanism produces the variable emission occurs similarly in all metal-poor 
red giants. 



The percentage of stars showing inflow and outflow in the Ha emission wingaj varies from 
cluster to cluster and appears to be related to cluster metallicity. In the metal-poor M92, about 
82% of s tellar spectra showing emission have an inflow signature (18% show outflow) and the study 
of M15 (IMeszaros et al.1 l2008l ) revealed that about 78% of stars with Ha emission displayed an 
inflow signature (22% outflow). These two clusters show similar behavior in their chromospheric 
dynamics. M13 has a more equal distribution of the dynamical signature with 55% of the Ha 
spectra indicating inflow (45% outflow). Since all the luminous stars are losing mass, it might 
be puzzling why the dominant emission signature in Ha is one of inflow. (And as discussed in 
the following section, the line cores generally indicate outflow.) We believe this relates to th e 
dynamics of the atmosphere. It appears likely that our targets are pulsating (jMavor et al.lll984! ). 
and comparison with the well-known pulsators, Cepheids, s hows that in flow signatures in Ha are 
accompanied by inflow in the metallic photospheric lines (jBaldryl 119971 ). and that inflow occurs 
for a larger fraction of the pulsation period in longer period Cepheids than in shorter period stars 



( Nardetto et al 



2006 



Petterson et al.ll2005l ). In this way, we understand the dominance of inflow 



emission signatures. The different proportions of inflow/outflow signatures may indicate that the 
pulsation period in M15 and M92 is generally longer than in M13. It is clear that variability is 
ubiquitous. Almost all stars brighter than V=12.5 are variables in M13, but only one variable red 
giant exhibits period ic photometric variations, and the remaining ones are semi-regular or irregular 
(|Kopacki etaD 120031 ) . 



The fraction of stars showing Ha emission increases with luminosity and decreasing effective 
temperature. Because of the separation of red giant branches in the CMD due to metallicity, the 
distribution of the emission with luminosity and effective temperatures differs between clusters. 
At the same luminosity, the metal rich M13 exhibits more stars with Ha emission than the metal 
poor M92 or M15, because the effective temperatures are lower in M13 for a constant luminosity. 



In M15, the stars K672 and K875 are clearly AGB objects, and the faintest RGB star showing emission in that 
study is B30. 

4 As measured by the ratio of B/R - the short wavelength to long wavelength peaks of the Ha emission wings. 
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However, at the same stellar effective temperature, M15 and M92 exhibit more stars showing 
emission than M13. This may appear as a metallicity effect, but it originates in the fact that both 
high luminosity and low effective temperature produce more Ha emission. 



6.2. The Ha Bisector Velocity 



Stars brighter than log (L/Lq) = 2.5 show a blue-shifted Ha core in both M13 and M92, 
and outflows become faster with increasing luminosity (Figure 9). This occurs in M15 as well 
(jMeszaros et al.ll2008l ). Thus, the luminosity at the onset of outflow, indicated by the Ha line core 
is independent of metallicity. The behavior of the bisector velocity on the RGB changes at the 
highest luminosities (Figure 14). Giants in M13, the most metal-rich cluster, show lower bisector 
velocities in the most luminous (and coolest stars). In fact, the velocities of the Ha core approach 
km s _1 with respect to the photosphere, signaling that the outward motions have decreased in 
the atmospheric region where the Ha line forms. Since the brightest stars in M13 have a lower T e ff 
than those in M15 and M92 (Fig. 14), we suspect that the decrease in the Ha bisector velocity 
results from the changing structure of the very extended atmosphere. The Ha wing asymmetry 
and Ca II K asymmetry preponderantly signal outflow in the most luminous stars. It is noteworthy 
that outflow begins at a luminosity, log (L/L®) ~ 2.5, and as the stars become more luminous, 
emission wings occur in the Ha profile in our sample at log (L/L®) ~ 2.75 . We conjecture that 
the onset of pulsation marked by the observed outward motion leads to a warmer chromosphere 
producing emission wings in H-a. 

Differences in the coreshift between AGB and RGB stars can be seen where these stars are 
distinct on the CMD. Stars on the AGB, between log (L/L®) = 2.0 and 2.7, showed slightly larger 
bisector velocities than RGB objects in both clusters, although the values are most extreme in M92. 
AGB stars tend to have lower values of log g and smaller escape velocities in the chromosphere 
as compared to RGB stars, which makes them more sensitive to mass loss driving mechanisms 
resulting in faster winds. We can speculate that there is more heating in the hotter AGB stars; it 
may be that a putative magnetic field is stronger after the stars have been through helium burning 
enabling enhanced wave motions, heating, and acceleration in the chromosphere as compared with 
RGB stars. The extremes of outflow velocity on the AGB tend to be larger in the metal-poor 
clusters, M15 and M92, than in M13 (see Figure 10 and 14). 

There is no evidence here that the outflow velocity is slower at low metallicity. This suggestion 



resul ted from observations of OH /IR star s in the low metallicity Magellanic Clouds ([Marshall et al 



2004 ). and modeling of dust driven winds (jHellingll2000l ; IWachter et al.ll2008l) such as those identified 



in O mega Centauri by Spitzer Space Telescope observations (jBover et al.l 120081 ; iMcDonald et al 



2009). In fact, we find just the opposite. M15 and M92 have genera lly higher velocities than stars in 
M13 (see Figure 14, lower panels). No evidence for a 'super-wind' (|Renzinilll98ll ; lBowen Willson 



1991 



) in the sense of an abrupt high velocity outflow is present in our spectra, altho ugh the largest 



mass- loss rates might be expected for stars with lower T e ff than found in this sample (jWachter et al 
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2002). Even in the dusty red giants in M15 detected with the Spitzer Space Telesco pe (Bover et al 



2006), the bisector velocities have similar values as red giants without an IR excess (jMeszaros et al 
2008). These similarities suggest that mass loss and dust production are not correlated, and the 
triggering of dust production may be an episodic phenomenon. 

Three stars in M13 exhibit large (>2 km s _1 ) bisector-velocity changes between observations. 
The star at the lowest luminosity limit, L719, shows a 4.7 km s _1 bisector velocity difference, which 
is clearly visible on the spectrum (see Fig. 3) as the Ha emission disappeared. K342 and K658 
showed 2.2 km s _1 and 4.9 km s _1 changes respectively, but these stars are faint and the error 
due to the low signal-to-noise of the spectra is comparable to the velocity change. In M92, only 
2 stars show a large coreshift in Ha: II-6 is a very faint star and this difference is comparable to 
the error of our measurements; IX-12 is an AGB star and shows similar vus values as stars of the 
same luminosity. However between our 2 obse rvations separated by 2 days, the coreshift changes 
for other stars are relatively small. In M15 (jMeszaros et al.1 120081 ) . it is the AGB stars in the 
log (L/Lq) = 2.3 — 3.0 luminosity region that show large bisector velocity changes (3—7 km s _1 ) 
over a time span of one year or more. 

On the RGB, the velocities in M15 indicate that systematic outflow (more negative than 
—2 km s _1 ) in the Ha core occurs at luminosities, log (L/Lq) > 2.5. The velocity of the outflow 
increases with luminosity and only the brightest stars show slightly smaller outflow (Figure 9 
and 14). In M92, which is also very metal-poor, stars brighter than log (L/Lq) = 2.5 showed 
bisector velocities more negative than —2 km s , and only the brightest star shows smaller outflow 
velocities. If there are differences in mass loss between M15 and M92, the shapes of the Ha line 
profiles do not reflect this. Thus t he Ha line by itself cannot he lp to decide if mass loss is the 
'second parameter' in M15 and M92. iMcDonald &: van Loonl (|2007l ) found no significant correlation 
between core asymmetry and luminosity, when they examined Ha and Ca II IR triplet spectra of 
47 red giant stars near the RGB tip in 6 globular clusters. Above a certain luminosity the bisector 
velocity of Ha becomes small and motions are difficult to detect in this region of the atmosphere, 
independent of cluster metallicity. Possibly another diagnostic such as He I A 10830 or ultraviolet 
lines, for med higher in the atmosphere needs to be examined. A He I A10830 absorption line was 



detected ( Smith et al.ll2004l ) in the star L687 in M13, with an extension to —30 km s , suggesting 



that when the helium line becomes detectable (apparently for T e ff ~ 4650 K), it can give an 
indication of a faster wind. This He I A10830 velocit y is comparable to t he hi gher values of the Ha 
bisector velocities that are found here in AGB stars. iPilachowski et al.l (| 19961 ) classifies L687 as an 
AGB star. 



6.3. The Ca II Emission and Velocity 



The Ca II K2 emission appears at lower luminosities than Ha emi ssion in both c l usters 
[log (L/L Q ) = 1.92 for M13 and log (L/L@) = 1.96 for M92]. In M15 (|Meszaros et al.l boOfih 
the Ca II K luminosity limit agrees with the Ha emission limit [log (L/Lq) = 2.36], but the low 
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signal-to-noise ratio of those observations did not allow us to determine the presence of Ca II H&K 
emission in fainter stars. The lower luminosity limit of Ca II K emission does not appear to be 
dependent on the cluster metallicity. 

The number of stars with Ca II emission in both M13 and M92 exceeds the number of stars 
showing Ha emission. Stars with Ha emission generally have Ca II K emission, but not all stars 
with Ca II K emission show Ha emission. This difference is not unexpected since the regions of 
formation of the Ca II K core and the Ha emission wings are separated in the atmosp here of a giant . 
Models suggest that Ca II K emission forms lower in the atmosphere than Ha wings (IDupredll 986). 
Addi t ionally Ha shows var iations in asymmetries over the span of a few days (jCacciari &: Freeman 



19831 ; iMeszaros et al.ll2008l ) which could contribute to the differences. Some stars in both clusters 



were observed twice at Ca II. Changes in Ca II K emission were observed in two stars in M13 and 
four stars in M92 where the line profile of Ca II K, either changed asymmetry or the emission 
strengthened or weakened, or both. 

The outflow velocities of the Ca II K3 reversal are generally larger th an the bisector velocities 
of the Ha line for the same stars (Figure 9). Similar behavior was found bv lZarro fc Rodgersl (|1983l ) 
in Population I giants and supergiants, and they concluded from the similarity of line profiles that 
Ca II K line is formed higher in the atmosphere and the increased outward velocity reflects a mass- 
conserving outflow. While models of the Sun sugge st that the C a II K3 feature forms in a higher 
atmospheric region th an the core of the Ha line (|Ayrettl Il998l ) , some chromospheric models of 
metal-deficient giants ( Dupree 1986 : Mauas et al. 20od ) locate the approximate depth of formation 
of the Ca II K3 feature below that of the Ha core. These models would suggest the opposite 
conclusion from Population I stars , that the flow is decelerating in the upper atmosphere. Yet 
another model ([Dupree et al.lll992l ) for the metal-deficient giant HD 6833 finds the contribution 
function of Ca II K3 to lie above that of the Ha core and hence signal an accelerating outflow. Some 
ambiguity may exist in the definition of the region of formation, and in addition it can extend over 
a substantial height in the chromosphere. In some cases, the He I A 10830 line, clearly formed above 



both Ca II K and Ha s hows even higher outflow velocities in metal-deficient stars (lDupree et al 



19921 : ISmith et al.ll2004l ). so the accelerating outflow models appear preferable. 



6.4. Comparison to NGC 2808 



Red giants in another metal rich cluster, NGC 2808, were studied by ICacciari et al.1 (120041 ) 
and can be compared to M13. However, this comparison may be somewha t compromised since 
NGC 2808 has an extreme case of peculiar horizontal branch morpholog y (|Lee et al.1 120051 ) and 
a split main sequence with potentially 3 populations ( Piotto et al.l 120071 ). making it one of the 
most persuasive clusters for the existenc e of possible multiple stell ar populations including a super 
helium-rich component jhee et al.1 bood : b'Antona & Caloil VMM ). These features of NGC 2808 
make it quite different from M13 - the metal rich cluster in our sample. M13 has a more chemically 
homogeneous population, although it is conjectured to consist of predominantly second generation 
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stars (ID'Antona & Caloil 12008 ) . The average metallicity of NGC 2808 ([Fe/H] 
than M13 ([Fe/H]=-1.54), by a factor of 2.5. 



T.15) is higher 



The detection threshold for Ha emission on the RGB in NGC 2808, of log (L/L@) = 2.5 is 
fainter by ~0.2 magnitudes than the limits for M13, M15, and M92 (see Table 11). Since the 
appearance of emission in the Ha line varies with time, these limits are comparable, one with 
another. However, the percentage of red giants exhibiting emission is less at 52% than we find for 
M13, M15, and M92 where the value is about 80%. The atmospheres of the NGC 2808 giants may 
be at lower temperatures since, for the same input energy, radiation losses are greater due to the 
increased abundance of metals than in metal-poor objects. 

Differences arise in the Ha outflow signature indicat ed by the emis s ion w ings in NGC 2808 
where an exceptionally low percentage (at 7%) is found by lCacciari et al.1 (J200J) in contrast to the 
45% of red giants showing outflow in M13 and 18—22% in the metal-poor clusters, M15 and M92. 



Cacciari et al.l (|2004l ) measured significant core shifts (< — 2 km s _1 ) of Ha in 7 stars of their 
sample of giants in NGC 2808. These stars are brighter than log (L/Lq) = 2.9. It is interesting 
that the 3 most luminous stars in their sample had core shifts of 1 km s~ 1 or less, similar to our 
results for M13 (Figure 9). M4, another cluster of similar metallicity as NGC 2808 also did not 
have coreshifts (more negative than —2 km s -1 ) eith er in Ha or Na D in any of ~10 stars that have 
luminosities log (L/L & ) > 3.3 (jKemp Bates! Il995l ). 

The luminosity limit in NGC 2808 (jCacciari et al. for Ca II K emission lines is log (L/Lq) <■ 

2.6 which is higher than the Ha limit in NGC 2808. This result is puzzling since Ca II K is found 
at lower luminosities than Ha emission in the other c l usters , M13, M15, and M92. The resolu- 
tion of the Calcium spectra studied by ICacciari et al.l (J200J) was the lowest of all their spectra 
at R=19600, and the signal-to-noise in the line core for the brightest stars was only 15. So it is 
possible that Calcium was not detected in the fainter targets. The limit for Ca II K in NGC 2808 
is ~0.2 magnitudes brighter than found in M15 which is a metal-poor cluster. The 2 metal poor 
clusters, M15 and M92 differ in the Ca II K limit by 0.4 magnitudes. At present, there is not 
sufficient evidence to claim that the B/R ratio of Ca II K emission varies night to night as the B/R 
ratio of Ha emission. The core shift of Ca II K in NGC 2808 is generally more negative than the 
value for Ha, similar to that found here for the most luminous stars (Figure 9). 

Since both the Ha emission and emission wing asymmetries are variable, it is difficult to draw 
firm conclusions abo ut systematic d i fferen ces between clusters. In the sample of red giants in 
NGC 2808 studied bv lCacciari et al.l (|2004l ). a lower fraction of stars was found with Ha emission 
and with outflow signatures in the emission wings than in the more metal-poor clusters (M13, M15, 
M92). However, the dynamical characteristics including the luminosity onset of outflow and wind 
speeds, appear indistinguishable among these clusters. 
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7. Conclusions 

Summarizing, we find the following: 

1. Hectochelle spectra of M13, M15, and M92 show Ha emission to occur on the red giant 
branch in stars with T e ff < 4500 K and log(L/L ) > 2.75. AGB stars exhibit Ha emission to 
lower luminosities. Ca II K emission extends to lower luminosities than Ha both on the RGB and 
AGB. 

2. Considering 3 clusters, spanning [Fe/H]=-1.54 (M13), to [Fe/H]=-2.3 (M15, M92), we 
find no systematic dependence of the presence of Ha or Ca II K emission from red giants on T e ff, 
L/L©, or cluster metallicity. 

3. Asymmetric Ha cores show that chromospheric material is flowing out from stars brighter 
than log (L/Lq) ~ 2.5 and the speed of the outflow increases with increasing stellar luminosity. 
This outflow may represent the onset of mass loss, and the luminosity at which outflow begins is 
similar for all metallicities. Ha velocities on the red giant branch are similar for all metallicities 
(but not for AGB stars, see below). 

4. Stars on the asymptotic giant branch near log (L/L©) ~ 2.0 — 2.7 show higher outflow 
velocities than RGB stars, and faster outflows are found in the metal— poor M15 and M92 than the 
metal— rich M13 objects. 

5. The sensitivity of Ha to mass motions decreases for T e ff < 4000 K causing the coolest 
giants in M13 to exhibit little or no outflow in this line. 

6. The Ca II K3 absorption features exhibit higher velocities than Ha suggesting accelerating 
outflows in the chromospheres. 

7. We find no differences in chromospheric signatures in the profiles or the presence of Ha and 
Ca II that can resolve the 'second-parameter' problem for the paired clusters, M15 and M92. 
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Fig. 1. — Color-magnitude diagram for all stars observed in M13 and M92. The solid line shows the 
fiducial curve of the RGB; the dashed line shows the fiducial curve of the AGB for both clusters 
taken from observations of ISandagej (|1970l ). The absolute magnitudes were calculated usi ng the 
appar ent distance modulus (m — M)y = 14.48 for M13 and (m — M)y = 14.64 for M92 from lHarris 
(Il996h . 
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Fig. 2. — Examples of the effect of subtracting an averaged spectrum (dashed line) from the observed 
spectra (solid lines) for stars in M13. The difference spectrum is shown by a dotted line below. L549 
is a star without any emission; the error of the subtracted spectrum (dotted line) is smaller than 
0.02 of the continuum level. In the case of L403, weak emission on the short wavelength is visible, 
however it is comparable to the noise of the observed spectrum and extends to the core of the line, 
so it was not identified as emission. L252 is an example of how the continuum normalization can 
shift the region near Ha making it hard to identify the emission. The blue emission in the spectrum 
of L250 might be missed by eye, but the subtraction method clearly shows the presence of emission. 
L465 and L72 are examples of emission that is clearly visible in the spectrum above the continuum 
level. 
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Fig. 3. — Normalized spectra of red giants in M13 which showed emission in Ha on 2006 March 
14. Stars marked with * were observed on 2006 May 10. The dashed line marks the bisector. The 
emission of one star, L719, disappeared between observations, and the spectrum is overlaid here 
using a dotted line. The spectra are arranged in order of decreasing brightness; the brightest is at 
the top left and the stars become fainter from left to right. The wavelength scale is corrected for 
heliocentric velocity. The radial velocity of M13 is —243.5 ± 0.2 km s _1 . 
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Fig. 4. — Normalized spectra of red giants in M92 which showed emission in Ha on 2006 May 7. Stars 
marked with * were observed on 2006 May 9. VII-10 is an example of an Ha profile without emission. 
For explanation please see Figure 3 caption. The radial velocity of M92 is —118.0 ± 0.2 km s -1 . 
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Fig. 5. — Color-magnitude diagrams for all M13 and M92 stars observed in 2006. Stars with Ha 
emission and with B<R (indicating outflow) are marked with filled circles; stars with B>R emission 
wings (suggests inflow) are denoted by filled squares. The solid line shows the fiducial curve of the 
RGB; dashed lines show the fiducial curve of the AGB for M13 and M92 from observations of 
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Fig. 6. — Continuum normalized spectra of a sample star (L96) in M13 showing Ha, RV31, and 
Ca II H&K spectra after all reductions. Upper spectrum is the observed one, lower spectrum is the 
model synthesis of a star (Coelho et al. 2005, using Kurucz models) with the highest amplitude of 
the cross-correlation function from the Ha region. The cross-correlation region used in the OB25 
filter is marked in the spectrum and chosen to avoid Ha. 
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March 14 (v ra a t i) compared to the same stars observed by ISoderberg et al.l (|1999l ) (v so d)- There 
is a slight offset (1.1 ± 0.5 km s -1 ) between all observations taken in 2006 and observations for 



the same stars from ISoderberg et al.l (jl999l ). Top right: Radial velocities measured in this paper 
on 2006 March 14 (v ra d,i) for the same stars observed on 2006 May 10 with the Ha filter (iv a d,2)- 
Lower left: Radial velocity measured with Hectochelle for the same stars observed on 2006 March 
14 (iVadl) compared to the observations with the RV31 filter on 2006 March 16 (v ra d,3)- Lower 
right: Radial velocities for the same stars measured with Hectochelle on 2006 March 14 (v ra d,i) 
compared to the observations with the RV31 filter on 2006 May 10 (v m <2,4)- The dashed line marks 
a 1:1 relation in all panels. The error of our measurements was generally smaller than the symbols 
used in the figure. The anomalous star in M13, L719, lies between the AGB and RGB, and the 
large velocity change may indicate binarity. 
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Fig. 8. — Top left: Radial velocities meas ured with the Ha filte r in this paper on 2006 May 7 {v ra d,i) 
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Fig. 9. — The velocity difference (ubis) between the top and the bottom of the bisector of Ha 
(x) and the coreshift of the Ca II K central reversal absorption (filled circle) as a function of 
luminosity. All Ca observations from different days are plotted together on the left side panels. 
Negative values indicate a blueshifted core (outward motion), positive values denote a red shifted 
core (inward motion). The error bars in figures on the left side were eliminated to display the 
differences between Ha and Ca II K. A predominant outward motion sets in near logL/L & « 2.5 
in both clusters and increases in velocity towards higher luminosity. The velocity of the Ca II K 
central reversal formed higher in the chromosphere than the Ha core, is generally larger than the 
bisector velocity of Ha at the same luminosity, indicating that the expansion velocity increases 
with height in the chromosphere. 
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Fig. 10. — Temperature-luminosity diagram for all stars observed in both clusters, where the size 
of the circle indicates the velocity of the Ha bisector asymmetry. Big circle: v^s < —4 km s _1 , 
medium circle: —4 km s _1 < vus < ~ 2 km s _1 , small circle: vu s > —2 km s _1 . Concentric circles 
indicate multiple observations of the same star. 
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Fig. 11. — Spectra of the brightest red giants in M13 which showed emission in Ca II K on 2006 
March 14. The spectra are smoothed by 3. The spectra are arranged in order of decreasing bright- 
ness; the brightest is at the top left and the stars become fainter from left to right for a single 
date. The object names marked by stars were observed on 2006 May 10. The wavelength scale is 
corrected for heliocentric velocity. The line marked ISM in the spectrum of L598 denotes absorption 
by the interstellar medium and can be recognized in the other spectra. 
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Fig. 12. — Spectra of the brightest red giants in M92 which showed emission in Ca II K on 2006 
May 7. The object names marked by stars were observed on 2006 May 9. Observations obtained on 
2006 May 7 have generally lower S /N due to bad sky conditions and resulted in a higher Ca II K 
core of XII-5, the faintest star in our sample. Additional explanation can be found in the caption 
of Figure 11. 
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Fig. 13. — Color-magnitude diagrams for all M13 and M92 stars observed. Stars with Ca II K 
emission and with B<R (indicating outflow) are marked with circles; stars with B>R emission 
wings (suggests inflow) are denoted by squares and stars with B~R are marked with triangles. The 
solid line shows the fiducial curve of the RGB; dashed lines show the fiducial curve of the AGB for 



M13 and M92 from observations of ISandagd (|1970l ). 
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Fig. 14. — Top left and right: The bisector velocity (v b i s ) of Ha for all observations in M13, M15 
and M92 as a function of luminosity and effective temperature. Negative values indicate outflow. 
Lower left and right: To fit the results with a linear function, luminosity and effective tempera- 
ture were divided into three different regions: the bottom of the RGB [log (L/L Q ) = 1.6 — 2.5, 
T e ff =4750 -5700 K], where v bis was close to zero km s _1 ; the RGB stars [log (L/L Q ) = 2.5-3.1, 
Tg//=4300— 4750 K], where vu s shows a significant increase; and the top of the RGB [log (L/L & ) = 
3.1 — 3.5, T e jj=3800— 4300 K], where v^s generally shows smaller values than in the middle of the 
RGB. The stars with high velocity near log (L/Lq) = 2.2 — 2.4 were omitted from the fit. Error 
bar of the top of the RGB fit for M13 is displayed in the lower left panel (for details see Section 
6.3). The errors of the fits span between ~ 1 and 2.5 km s _1 for each fit. 
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Table 1. Photometric Data of Observed Cluster Members in M13 



ID No. a 


RA(2000) b 


Dec(2000) b 


B 


V 


,1 


H 


K 


B-V 


V-K 


Obs. c 


111-65 d 


16 41 39.091 


+36 23 51.40 


16.04 


15.22 


13.690 


13.220 


13.071 


0.82 


2.15 


1,3,5 


K188 c 


16 40 42.982 


+36 27 41.88 


14.44 


13.39 


11.410 


10.784 


10.704 


1.05 


2.69 


1,3,5 


K210 ° 


16 40 56.378 


+36 22 18.51 


15.21 


14.33 


12.580 


12.028 


11.966 


0.88 


2.36 


1,3,5 


K220 ° 


16 41 02.608 


+36 26 15.81 


15.74 


14.98 


13.353 


12.899 


12.833 


0.76 


2.15 


1,2,4,5,6 


K223 ° 


16 41 05.075 


+36 28 20.85 


15.56 


14.71 


13.093 


12.565 


12.494 


0.85 


2.22 


2,4,6 



3 Ludcndor | h905h is the identification for the majority of the stars denoted by L. 
b 2MASS coordinates llSkrutskie et alj|2006ri 

Observations: 1: 2006 March 14 (OB25), 2: 2006 May 10 (OB25), 3: 2006 March 16 (RV31), 4: 2006 May 10 (RV31), 
5: 2006 March 16 (Ca41), 6: 2006 May 10 (Ca41). 

' Urd dl955h 

jKadlal jl966h 

Note. — The visual photome try is taken from lCudworth &: Monet] lll979l) , J,H,K photometry is taken from the 2MASS 
Catalog l lSkrutskie et alj|2006h . This table is available in its entirety in a machine-readable form in the on-line journal. 
A portion is shown here for guidance regarding its form and content. 



Table 2. Photometric Data of Observed Cluster Members in M92 



ID No. a 


RA(2000) b 


Dec(2000) b 


B 


V 


,1 


H 


K 


B-V 


V-K 


Obs. c 


1-14 


17 17 28.77 


43 10 02.8 


15.47 


14.74 


13.155 


12.644 


12.592 


0.73 


2.148 


1,2,3,4,5 


1-40 


17 17 22.68 


43 08 50.5 


15.51 


14.78 


13.258 


12.777 


12.640 


0.73 


2.140 


1,3,4 


1-67 


17 17 21.24 


43 08 27.0 


14.24 


13.32 


11.406 


10.870 


10.766 


0.92 


2.554 


1,3,4 


1-68 


17 17 21.73 


43 08 15.8 


15.36 


14.61 


13.243 


12.825 


12.661 


0.75 


1.949 


2,5 


II-6 


17 17 50.37 


43 13 46.0 


15.89 


15.14 


13.541 


13.002 


12.992 


0.75 


2.148 


1,2,3,4,5 



HSandage fc Walkerl jl966h 

b 2MASS coordinates l|Skrutskie et al.ll2006ri 

Observations: 1: 2006 May 7 (OB25), 2: 2006 May 9 (OB25), 3: 2006 May 7 (RV31), 4: 2006 May 8 (Ca41), 5: 
2006 May 9 (Ca41). 



Note. — The visual photom etry is taken from ICudwortbl lll976h , J ,H,K photometry is taken from the 2MASS 
Catalog llSkrutskie et al.ll2006t) . This table is available in its entirety in a machine-readable form in the on-line 
journal. A portion is shown here for guidance regarding its form and content. 
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Table 3. Hectochelle Observations of M13 and M92 



Date 


Total exp. 


Wavelength 


Filter Name 


Number of 


(UT) 


(s) 


(A) 




Observed Stars 


2006 March 14 (M13, Field 1) 


3 x 2400 


6475-6630 


OB25 


70 


2006 March 16 (M13, Field 1) 


3 x 2400 


3910-3990 


Ca41 


70 


2006 March 16 (M13, Field 1) 


1 x 2400 


5150-5300 


RV31 


65 


2006 May 10 (M13, Field 2) 


3 x 2400 


6475-6630 


OB25 


70 


2006 May 10 (M13, Field 2) 


3 x 2400 


3910-3990 


Ca41 


63 


2006 May 10 (M13, Field 2) 


1 x 2400 


5150-5300 


RV31 


65 


2006 May 7 (M92, Field 1) 


3 x 2400 


6475-6630 


OB25 


42 


2006 May 7 (M92, Field 1) 


3 x 1800 


5150-5300 


RV31 


40 


2006 May 8 (M92, Field 1) 


3 x 2400 


3910-3990 


Ca41 


41 


2006 May 9 (M92, Field 2) 


3 x 1800 


6475-6630 


OB25 


36 


2006 May 9 (M92, Field 2) 


3 x 2400 


3910-3990 


Ca41 


36 



-36- 



Table 4. B/R ratio of Ha Line for Stars with Emission Wings 



ID No. 


M13 
B/R 
2006 March 14 


M92 

B/R B/R B/R 
2006 May 10 ID No. 2006 May 7 2006 May 9 


L70 


< 1 


11-53 > 1 > 1 


L72 


< 1 


> 1 111-65 > 1 


L96 


< 1 


IV-94 > 1 


L158 


> 1 


VII-18 < 1 < 1 


L169 


> 1 


VII-122 > 1 


L199 


< 1 


IX-12 >1 >1 


L250 


> 1 


X-49 > 1 


L252 


> 1 


XII-8 > 1 


L316 


> 1 




L345 




> 1 


L414 


< 1 




L465 


> 1 




L598 


< 1 




L719 


< 1 


no emission 


L745 


> 1 




L835 


< 1 




L954 


< 1 




L973 


> 1 




L1073 


> 1 




Note. — The parameter B/R is the intensity ratio of Blue (short wavelength) and Red 
(long wavelength) emission peaks. The symbol ■ ■ ■ indicates the star was not observed. 
If B/R ratio is > 1 the emission wings indicate inflow, if B/R ratio is < 1 the emission 
wings indicate outflow. 
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Table 5. Physical Parameters of Cluster Members in M13 



ID No. 


M v 


(B - V) 


(V-K) 


P a 


T eff 

(K) 


log L/Lq 


R/R Q 


111-65 


+0.74 


0.80 


2.095 


99 


5010 


1.750 


9.7 


K188 


-1.09 


1.03 


2.635 


99 


4470 


2.575 


31.5 


K210 


-0.15 


0.86 


2.305 


99 


4790 


2.138 


16.6 


K220 


+0.50 


0.74 


2.095 


99 


5010 


1.846 


10.8 


K223 


+0.23 


0.83 


2.165 


99 


4900 


1.964 


13.0 



"■M embership probability from proper motion observations dCudworth fc Monetl 

Eizi). 



Note. — This table is available in its entirety in a machine-readable form in the 
on-line journal. A portion is shown here for guidance regarding its form and content. 



Table 6. Physical Parameters of Cluster Members in M92 



ID No. M v (B-V)o (V-K) P a T e// log L/L e R/R Q 

(K) 



1-14 


+0.10 


0.71 


2.093 


99 


4980 


2.006 


13.5 


1-40 


+0.14 


0.71 


2.085 


99 


4990 


1.989 


13.2 


1-67 


-1.32 


0.90 


2.499 


99 


4570 


2.640 


33.4 


1-68 


-0.03 


0.73 


1.894 


99 


5240 


2.033 


12.6 


II-6 


+0.50 


0.73 


2.093 


99 


4980 


1.846 


11.3 



"Membership probability from proper motion observations l|Cudworthlll97l 



Note. — This table is available in its entirety in a machine-readable form in the 
on-line journal. A portion is shown here for guidance regarding its form and content. 
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Table 7. Radial and Ha Bisector Velocity of Observed Stars in M13 



ID No. 


V ra d,l a 


*W,2 a 


V r ad,3 a 


v rad,4 




Vbis,2 b 




(km s" 1 ) 


(km s" 1 ) 


(km s" 1 ) 


(km s- 1 ) 


(km s" 1 ) 


(km s- 1 ) 


111-65 


-246.0 ± 0.4 




-245.9 ± 0.2 




+0.5 ± 2.5 




K188 


-244.9 ± 0.3 




-245.2 ± 0.2 




-1.4 ± 1.0 




K210 


-250.8 ± 0.2 




-249.9 ± 0.3 




-0.0 ± 0.6 




K220 


-244.0 ± 0.3 


-243.6 ± 0.2 




-242.9 ± 0.2 


+0.4 ± 1.4 


+1.5 ± 0.7 


K223 




-242.6 ± 0.2 




-242.0 ± 0.2 




+0.2 ± 0.9 



a Observations: 1: 2006 March 14 (Ha), 2: 2006 May 10 (Ha), 3: 2006 March 16 (RV31), 4: 2006 May 
10 (RV31). 

Observations: 1: 2006 March 14; 2: 2006 May 10. 

Note. — This table is available in its entirety in a machine-readable form in the on-line journal. A 
portion is shown here for guidance regarding its form and content. 



Table 8. Radial and Ha Bisector Velocity of Observed Stars in M92 



ID No. 


V r ad,l a 


Vrad, 2 3 


Vrad,3 ^ 




Vbis,2 b 




(km s" 1 ) 


(km s" 1 ) 


(km s" 1 ) 


(km s" 1 ) 


(km s" 1 ) 


1-14 


-124.2 ± 0.3 


-124.0 ± 0.2 


-123.9 ± 0.2 


-0.5 ± 1.2 


+0.6 ± 0.5 


1-40 


-125.2 ± 0.2 




-124.6 ± 0.2 


+0.5 ± 0.5 




1-67 


-120.9 ± 0.2 




-120.2 ± 0.2 


-1.8 ± 0.4 




1-68 




-129.7 ± 0.2 






+0.0 ± 0.6 


II-6 


-122.6 ± 0.3 


-123.1 ± 0.3 


-121.8 ± 0.2 


+2.5 ± 0.7 


-1.8 ± 0.8 



a Observations: 1: 2006 May 7 (OB25), 2: 2006 May 9 (OB25), 3: 2006 May 7 (RV31) 
a Observations: 1: 2006 March 14; 2: 2006 May 10. 



Note. — This table is available in its entirety in a machine-readable form in the on-line 
journal. A portion is shown here for guidance regarding its form and content. 
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Table 9. B/R ratio of Ca II K Line for Stars Showing Emission in M13 and M92 



M13 M92 
B/R B/R B/R B/R 

ID No. 2006 March 16 2006 May 10 ID No. 2006 May 8 2006 May 9 



K228 


> 1 




11-53 


> 1 


S3 1 


K422 


S3 1 


S3 1 


111-65 


> 1 




K656 


> 1 




VII-10 


S3 1 


S3 1 


L18 


> 1 




VII-18 


< 1 


< 1 


L26 


> 1 


S3 1 


VII-39 


< 1 




L70 


< 1 




VII-122 


> 1 




L72 


< 1 




X-49 


> 1 




L95 




> 1 


XI- 14 


S3 1 




L96 


< 1 




XI-80 


> 1 




L109 


> 1 




XII-5 


> 1 


S3 1 


L140 


> 1 




XII-8 


> 1 




L158 


S3 1 




XII-34 


S3 1 




L169 


< 1 










L199 


< 1 










L250 


S3 1 










L252 


> 1 










L316 


> 1 










L345 




> 1 








L384 




< 1 








L403 


> 1 










L414 


< 1 










L423 


S3 1 










L465 


> 1 










L549 


S3 1 










L598 


S3 1 










L745 


S3 1 










L773 


> 1 










L835 


S3 1 










L863 


> 1 










L954 


< 1 










L973 


< 1 










L1023 


> 1 










L1043 


S3 1 










L1073 


S3 1 











Note. — The parameter B/R is the intensity ratio of Blue (short wavelength) and Red 
(long wavelength) emission peaks. 
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Table 10. Relative Radial Velocity of Ca II K central absorption. 



M13 M92 



ID No. 


V r el 


ID No. 


■"re.1 




(km s- 1 ) 




(km s- 1 ) 


K228 


+0.2 ± 0.8 


11-53 


-7.2 ± 0.9 


K422 


-3.2 ± 0.9 


111-65 


-10.9 ± 0.9 


K656 


-2.2 ± 0.7 


VII-18 


-13.0 ± 0.4 


L18 


-0.3 ± 0.8 


VII-122 


-6.6 ± 0.8 


L26 


-3.5 ± 0.9 


X-49 


-7.6 ± 1.0 


L70 


-10.2 ± 0.5 


XII-8 


-5.4 ± 0.8 


L72 


-10.0 ± 0.9 






L96 


-14.7 ± 0.9 






L109 


-7.8 ± 1.0 






L140 


—0.6 i 0.6 






L158 


-7.1 ± 0.9 






L169 


-6.4 ± 0.9 






L199 


-10.8 ± 0.8 






L250 


-12.8 ± 1.0 






L252 


-1.4 ± 0.8 






L316 


-6.1 ± 0.7 






L403 


-1.6 ± 0.9 






L414 


-12.0 ± 1.1 






L423 


-5.2 ± 1.1 






L465 


-6.3 ± 0.9 






L549 


-3.7 ± 1.2 






L598 


-7.2 ± 0.2 






L745 


-8.7 ± 0.8 






L773 


-3.9 ± 1.1 






L835 


-11.5 ± 0.9 






L863 


-4.0 ± 1.1 






L954 


-9.8 ± 1.1 






L973 


-11.0 ± 0.6 






L1023 


-3.4 ± 1.1 






L1043 


-0.8 ± 0.9 






L1073 


-9.4 ± 1.0 







Note. — The table does not contain stars where the central absorp- 
tion was not visible in the spectrum. 
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Table 11. Characteristics of Emission in Four Clust 


ers 






Cluster 


[Fc/H] a 


log(L/L QtHaA ) h 


log(L/LQ :Ha ,2) c 


log(L/L e>Ca k) d 


No. ° 


Pi ! 


P 2 g 


M13 


-1.54 


1.95 


2.79 


1.92 


123 


45 


78 


M15 


-2.26 


2.36 


2.78 


2.36 


110 


22.5 


83 


M92 


-2.28 


2.05 


2.74 


1.96 


64 


18 


78 


NGC 2808 h 


-1.15 


2.5 


2.5 


2.60 


137 


7 


52 



^Harris! [|l996ri . 

b The luminosity limit of all stars showing emission in Ha. 
c The luminosity limit of only RGB stars showing emission in Ha. 
d The luminosity limit of all stars showing emission in Ca II K. 
c Number of stars observed. 

f The percentage of stars from all observations showing outflow in Ha emission wing asymmetry. 

g The percentage of stars from all observations showing Ha emission above log(L/ Lq Ha,2)- 

h Parameters of NGC 2808 for the RGB stars were measured by ICacciari et al. No AGB stars in 

that sample are well-separated on the CMD. Thus the RGB limit marks the faint limit for all stars in the 
cluster. 



